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SUMMARY 


The  maximum  likelihood  equations  used  to  determine  estimates 
of  the  mean  and  standard  deviation  for  the  Bruceton  method  of 
statistical  analysis  were  solved  numerically.  A  digital  computer 
algorithm  was  developed  to  reduce  the  time  associated  with  the 
numerical  calculations.  The  method  of  successive  approximations 
(Newton- Raphson  method)  was  suitable.  Test  data  was  evaluated  to 
determine  the  relative  difference  between  values  of  the  means  and 
standard  deviation  as  obtained  from  the  numerical  solution  of  the 
maximum  likelihood  equations  and  as  obtained  from  formulae 
approximating  these  equations.  Results  of  the  analysis  showed 
that  the  differences  between  means  were  slight  while  those  between 
standard  deviations  were  more  pronounced  -  Calculated  high  and 
low  probability  levels,  because  of  their  dependence  on  the  mean  and 
standard  deviation,  differed  widely  from  one  method  of  calculation 
to  the  other.  The  analysis  also  showed  that  the  most  adverse  effects 
occurred  for  the  smallest  samples.  Because  the  maximum  likeli¬ 
hood  equations  yield  better  estimates  of  the  mean  and  standard 
deviation  than  can  be  obtained  from  equations  approximating  them  -- 
especially  for  small  samples  --  solution  of  the  maximum  likelihood 
equations  by  the  Newton- Raphson  method  is  advocated  for  determining 
the  best  estimates  of  the  population  mean  and  standard  deviation. 
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CONCLUSIONS 


& 


t 


The  maximum  likelihood  equations  used  to  estimate  the  mean 
and  standard  deviation  in  conjunction  with  a  Bruceton  analysis  are 
readily  solvable  by  the  Newton- Raphson  method  of  successive 
approximations  with  the  aid  of  the  digital  computer  codes  of 
Appendix  A  and  B. 

The  code  of  Appendix  B  appears  to  be  superior  based  on  higher 
reliability  and  safety  criteria. 
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RECOMMENDATION 


Because  of  its  apparent  superiority,  the  code  of  Appendix  B 
should  be  used  to  calculate  the  maximum  likelihood  estimates  of 
the  mean  and  standard  deviation  for  samples  analyzed  by  the 
Bruceton  method. 
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BACKGROUND 


The  Bruceton  method  of  statistical  analysis  is  used  to  determine 
the  probability  of  explosion  or  non-explosion  for  a  given  explosive 
item  as  a  function  of  input  stimuli.  The  procedures  used  for  calcu¬ 
lating  the  statistical  parameters  {mean,  standard  deviation, 
percent  points,  etc.)  are  based  on  approximations  to  the  maximum 
likelihood  equations  which  will  be  discussed  at  greater  length  in  the 
Study  section.  However,  the  validity  of  the  approximations  fall  off 
rapidly  with  decrease  in  sample  sizes  less  than  100.  Consequently, 
it  is  desirable  to  use  the  maximum  likelihood  equations  for  small 
sample  analysis. 

Before  the  advent  of  the  digital  computer,  exact  solutions  of 
the  maximum  likelihood  equations  required  complex,  wholly 
impractical  computational  schemes.  As  a  result  of  the  difficulties 
formerly  encountered  in  applying  the  likelihood  theory,  approximate 
methods  were  always  used  to  estimate  the  maximum  likelihood 
quantities.  In  the  case  of  small  sample  analysis,  blatant  deviations 
frequently  occurred  between  sample  and  population  statistical 
parameters.  The  Franklin  Institute  of  Philadelphia,  Pennsylvania 
investigated  the  affect  of  approximate  methods  on  the  data  obtained 
from  small  samples  and  their  findings  indicate  that  the  expected 
deterioration  of  the  approximate  formulae  as  applied  to  small 
samples,  does  in  fact  occur. 

It  must  be  emphasized  that  in  any  statistical  analysis  large 
samples  are  preferable.  But  if  small  samples  must  be  used  in 
a  Bruceton  test,  the  use  of  maximum  likelihood  techniques  become 
in  a  realistic  sense,  mandatory.  Furthermore,  the  maximum 
likelihood  vch'.es  associated  with  small  samples  should  be  treated 
with  care,  since  small  samples  tend  to  yield  poor  statistical 
results.  This  is  especially  true  for  this  particular  problem. 
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STUDY 


A  method  of  statistical  analysis  frequently  utilized  in  evaluating 
explosive  devices  is  known  as  the  Bruceton  Staircase  Sensitivity 
Test.  The  method  was  developed  by  the  Explosives  Research 
Laboratory  cf  Bruceton,  Pennsylvania.  It  is  used  to  ascertain  the 
safe  and  sure-fire  energy  levels  associated  with  a  given  explosive 
item.  In  performing  the  experimental  procedure  the  tester  is  first 
required  to  home  in  on  a  50%  firing  level  (mean)  and  on  an 
approximate  standard  deviation  (interval  between  consecutive  energy 
levels).  This  can  be  accomplished  easily  if  the  past  history  of  the 
explosive  under  investigation  is  available.  Otherwise  trial  and 
error  techniques  are  adequate.  Once  having  chosen  the  50%  firing 
level  and  the  step  interval  (based  on  an  approximation  of  the  standard 
deviation)  the  tester  may  then  proceed  as  follows: 

An  attempt  is  made  to  actuate  the  explosive  using  the 
50%  energy  estimate.  If  successful  a  second  specimen  is  tested 
at  an  energy  level  one  step  interval  below  the  50%  value.  On  the 
other  hand,  if  the  specimen  failed  to  detonate,  the  second  explosive 
item  is  tested  at  an  energy  level  one  step  interval  above  the  50% 
level.  Thus  denoting  any  test  level  as  H|  and  the  step  interval  as 
d,  then  -  d,  for  the  case  where  the  i'  th  test  produced 

an  explosion  and  +  d  if  the  i1  th  test  produced  a  non¬ 

explosion. 

In  this  way,  a  sequence  of  explosions  (xs)  and  non-explosions 
(0s)  will  be  generated.  Figure  1  shows  a  typical  sequence.  It  may 
be  noted  that  the  bulk  of  the  testing  occurs  in  the  neighborhood  of 
the  mean  energy  level.  Thus  a  minimal  number  of  tests  will  permit 
a  good  estimation  of  the  mean.  However,  because  of  the 
approximate  nature  of  the  methods  of  calculation  used  to  estimate 
the  mean  and  standard  deviation,  samples  of  at  least  100  specimens 
are  required.  Also,  sensitivity  data  does  not  consist  of  exact 
measurements.  Consequently  the  information  contained  in  them  is 
reduced  by  about  one-half  that  which  would  be  contained  in 
measured  data  of  a  comparable  nature. 
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To  reduce  the  bias  associated  with  the  Bruceton  method, 
large  sample  sizes  (at  least  100  specimens)  should  be  used. 

For  very  expensive  explosive  devices,  large  sample  sizes  are 
prohibitive.  Consequently,  tests  are  conducted  using  smaller 
sample  sizes  often  50  specimens  or  less.  But  it  is  known  that 
for  sample  sizes  below  50,  the  approximations  of  the  mean  and 
standard  deviation  become  rather  poor.  Therefore,  more 
elaborate  methods  using  the  maximum  likelihood  equations  for 
estimating  the  mean  and  standard  deviation  must  be  used. 

However,  these  equations  cannot  be  solved  in  closed  form  by 
any  direct  method  (Equation  47,  p.  38  of  Reference  4).  Of 
course,  numerically  the  equations  can  be  solved  and  the  maximum 
likelihood  estimates  of  the  mean  and  standard  deviation  determined. 

Owing  to  the  time  involved  in  the  numerical  solution  of  the 
maximum  likelihood  equations,  no  attempt  was  made  tc  utilize 
them  until  the  advent  of  the  digital  computer.  However,  once 
the  computer  algorithm  developed  in  this  report  was  established 
it  became  a  rather  simple  task  to  demonstrate  the  effect  of 
sample  size  on  the  calculated  mean  and  standard  deviation 
(: specially  since  it  is  this  parameter  which  strongly  affects  the 
outcome  of  the  extrapolation  calculations.  That  is,  the  validity 
of  the  approximation  formulae  for  estimating  high  and  low  percent 
points  is  doubtful  --  especially  if  small  sample  sizes  are  used. 
Therefore,  the  purpose  of  this  report  is  to  present  a  numerical 
procedure  for  solving  the  maximum  likelihood  equations  and 
suitable  for  application  to  a  digital  computer.  The  effect  of 
sample  size  on  the  results  of  both  the  approximate  and  the 
numerical  techniques  will  be  emphasized. 

To  determine  the  maximum  likelihood  equations  for  estimating 
the  mean  and  standard  deviation,  of  a  normally  distributed  random 
variable  x  with  mean  A  and  variance  S  ,  a  random  sample  of  size 
N  is  drawn  from  the  population  and  tested  in  the  manner  described 
so  that  a  sequence  of  ones  and  zeros  is  obtained.  If: 

Nj  =  total  number  of  ones 
Nq  =  total  number  of  zeros 
H^  =  energy  at  test  level  i 
=  number  of  ones  at 
M.^  =  number  of  zeros  at 
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then 

Ni  =  iN'i 

N0=  iM; 

and  N  -  NQ  +  Nj 

Further,  it  may  be  assumed  Nj  <  Ng,  then  the  probability  that  an  X 
will  be  a  one  at  the  level  is: 


and  the  probability  that  X  will  be  a  zero  is:  Qj  =  1  -  , 


Then,  the  probability  of  generating  a  sequence  of  Nj  ones  and  Ng 
zeros  is  obtained  from: 


P..77W* 

where  K  is  a  function  of  the  Nj  and  Mj,  This  probability  is  based 
on  generating  a  sequence  of  N  independent  outcomes  (ones  and 
zeros)  from  a  multi-nomial  distribution.  Although  this  assumption 
is  not  strictly  applicable  to  the  current  problem,  only  slight 
adjustments  need  be  made  to  the  general  theory  to  apply  it  in  the 
asymptotic  case.  This  is  denoted  by  the  likelihood: 


LOG  j 


and  differentiating  with  respect  to  A  and  S  maximized  likelihood 
equations  are  obtained.  Before  proceeding  further  the  following 
derivatives  are  defined: 


l 

S 


Zi 


9  Pi  Hj  -  A 

ds  "  s2 


Zi 
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d  Qj  i 

d  A  =  S  ^ 

and  P  Q,  Hj  -  A 

~dT  "  IT  Zl 

where  Z.  =  _ L_  ^  ^  '  A>2 

1  e  t?" 


Then,  the  maximum  likelihood  equations  for  estimating  the  mean 
and  standard  deviation  are: 


2A 


(i) 


and 


3h 

as 


?  I  -  N 


-  A)Zj_  !  (Hi-A)Z 


Ps 


-  -0  (2) 


These  equations  cannot  be  solved  by  direc*  ethods.  However,  if 
good  estimates  of  the  mean  and  standard  iation  are  available  a 
method  of  successive  approximations  (T  ..on-Raphson  method)  is 
applicable.  The  method  ;  1  to:  lows:  Let  the  maxi¬ 

mum  likelihood  Equation  \  ’./  c*- 


Fj  (S,  A)  =  0  (3) 

F2  (S,  A)  =  0  (4) 


Also,  let  Ao,  So  be  approximate  solutions  of  Equation  (3)  and  (4) 
such  that  the  true  roots  of  Equation  (3)  and  (4),  A  and  S  are 
related  to  Aq  and  Sq  by: 


S  =  S0  +  H 
A  =  AQ  +  K 
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where  H  and  K  are  the  corrections  for  Sq  and  Aq,  respectively. 
Rewriting  Equation  (3)  and  (4): 


Fj  (SQ  +  H,  Aq  +  K)  =  0  l 

f2  (S0  +  K,  Aq  +  K)  =  0  I 

Equation  (5)  and  (6)  may  be  expanded  by  Taylor's  theorem  for  a 
function  of  two  variables  giving: 


Fx  (S0  +  H,  A0  +  K)  =  0  =  Fj 


Terms  in  Higher  Powers  of  H,  K 


f2  (S0  +  h,  a0  +  k)  =  o  =  f2 


(S0,  Aq)  +  Hp  FZ)  +  Kid  FZ)  + 

*3 S / 0  U  A/o 


Terms  in  Higher  Powers  of  H,  K 

Since  the  values  of  H  and  K  are  small,  their  squares,  products, 
and  higher  powers  are  neglected  so  that: 


F1  <so 


• v  *  *  iPl  ■ ° 


«  •  »(|^)  *  ■  • 


Using  Cramer's  rule  (assuming  D  4  0)  solution  for  H  and  K 
are  obtained  from: 
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and 


DET 


lw\ 

(§~s)  o  -  Fz  (s  °>  ^*) 
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(10) 


Using  the  values  of  H  and  K  obtained  above  new  approximations  are 
gotten  for  S  and  A  as: 


S1  -  so  +  Ho 


A1  =  Ao  +  K0 


If  the  process  is  repeated  using  and  new  corrections,  Hj 
and  Kp  are  obtained.  Continuing  in  this  manner,  a  predetermined 

degree  of  convergence  may  be  achieved.  The  subscripts  0,  1 . . 

refer  to  the  values  A,  S,  K,  and  H  to  be  used.  That  is  (  ^  ^  1  ) 

3F)  ( 0 

means  that  ~ -  is  to  be  evaluated  for  S  =  Sq,  A  =  Aq,  and  so  on. 

d  s 

Also,  examination  of  Equation  (7),  (8),  (9)  and  (10)  shows  that  for 
large  D  and  correspondingly  small  H  and  K,  the  assumptions  made 
are  applicable  and  convergence  should  be  expected. 

In  accordance  with  this  method,  Equation  (1)  and  (2)  must  be 
differentiated  partially  with  respect  to  S  and  A.  However,  before 
proceeding  the  following  quantities  are  needed: 


Then  for  Equation  (1): 


2£, 

?s 


S 

/ 

P2  / 
5S 


I*  £&*(&)- ’■££>*. ' 


and 


£?£' 

<P/9 


^  W  -  p4^  *  'I 


i4 


i 


r 

■Sir  " 


Repeating  the  preceding  for  Equation  (2)  the  following  results  are 
obtained: 
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Using  Equation  (9)  and  (10)  repeatedly,  corrections  (H  and  K)  are 
calculated  for  S  and  A  until  the  desired  degree  of  convergence  is 
achieved.  Of  course,  owing  to  the  lengthy  calculations  required 
by  this  method,  a  digital  computer  program  for  working  out  the 
details  was  essential.  Otherwise,  the  method  would  prove 
impractical. 
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Two  similar  FORTRAN  codes  for  calculating  the  results  of  a 
Bruceton  analysis  are  in  Appendix  A  and  B.  The  two  programs 
are  basically  the  same.  A  comparison  of  them  follows. 

The  main  program  entitled  BRUCE  in  both  codes  is  used  to 
read  in  information  pertaining  to  reliability  and  confidence  levels. 
Then  subroutine  FOREST  is  called.  (This  subroutine  is  a 
modification  of  a  FORTRAN  code  written  by  Forrest  L.  McMains, 
(Reference  3).) 

Subroutine  FOREST  is  used  to  read  in  the  basic  data  generated 
in  the  Bruceton  analysis  and  then  calculate  estimates  of  the  mean 
and  standard  deviation  from  approximations  of  the  maximum  likeli¬ 
hood  equations.  However,  the  estimates  returned  to  BRUCE  differ 
for  the  two  codes.  Using. the  code  of  Appendix  A,  the  estimates  of 
A  and  S  are  calculated  for  both  the  number  of  fails  and  fires  --  the 
one  obtained  for  the  smaller  of  the  preceding  quantities  being  returned. 
Actually,  the  choice  is  made  in  BRUCE  after  returning  from  FOREST. 
A  subroutine  called  SWAP  is  called  when  the  number  of  fails  is  less 
than  the  number  of  fires.  In  this  event  the  number  of  fails  will  be 
retained  for  further  calculation.  If  the  number  of  fires  is  less  than 
or  equal  to  the  number  of  fails  no  transfer  to  subroutine  SWAP  will 
take  place.  Thus  all  calculations  using  the  code  in  Appendix  A 
depend  on  the  smaller  of  the  two  quantities,  number  of  fails  and 
number  of  fires.  This  approach  is  advocated  in  Reference  4. 

Using  the  code  of  Appendix  B,  the  number  of  fires  and  fails 
are  equally  weighted  regardless  of  which  is  smaller.  The  basis 
for  this  approach  is  described  in  Reference  3.  Both  the  number 
of  fires  and  fails  are  used  throughout  the  program.  Subroutine 
FOREST  computes  approximations  to  the  mean  and  standard 
deviation  by  first  calculating  estimates  of  these  quantities  based  on 
the  number  of  fails  and  fires  alone  and  then  averaging  them.  These 
averages  are  weighted  so  that  the  quantities  calculated  using  the 
smaller  value  between  the  number  of  fires  or  fails  is  more  heavily 
favored. 

Thus  using  either  code  approximate  estimates  of  A  and  S  are 
obtained  using  FOREST.  After  writing  out  the  information  read 
in  plus  the  starting  values  of  S  and  A  just  chosen,  a  subroutine 
entitled  PROB  is  called.  This  subroutine  calculates  the  probabili¬ 
ties  P^  and  0^  associated  with  each  energy  level  H^,  The  procedure 
followed  in  this  subroutine  is  a  modification  of  an  existing  program 
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used  for  calculating  the  probabilities  of  normally  distributed 
random  variables  (Reference  2).  The  program  which  was 
modified  is  used  to  calculate  probabilities  from  the  normal 
probability  integral  defined  as: 


The  integral  gives  the  probability  that  the  normalized  random 
variable  t  will  be  between  -t  and  t  on  the  normal  probability 
curve  shown  in  Figure  2.  That  is,  the  probability  P(x)  is  equal 
to  the  area  under  the  curve  from  -t  to  t.  For  small  t  the  integral 
above  can  be  evaluated  from  a  power  series: 


Fort  >  2.90  a  more  efficient  method  of  evaluating  P(x)  is  achieved 
by  using  the  asymptotic  series  defined  as: 


P(x)  = 


w» 


itt 


/•2  /'3  'S’ 

T7" 


- 


For  small  t,  (t  <  1.20),  the  economy  of  calculating  the  normal 
probabilities  is  enhanced  if  use  is  made  of  telescoping  the  power 
series.  For  all  terms  beyond  and  including  t®,  use  is  made  of 
the  telescoping  effect  produced  by: 


t8  =  —  (256t2  -  160t4  4-  32t2  -  1  +  T8) 

128  v 

which  is  deduced  from  Chebyshev  polynomials  (Reference  2). 

For  the  range  1.2  <  t  <  2.9  application  of  the  original  power 
series  is  satisfactory. 
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FIGURE  2 


For  each  range,  t  <  1.20,  t  >  2.90,  and  1,20  <  t  <  2.90, 
probabilities  were  calculated  and  compared  with  those  recorded 
in  statistical  tables.  Because  of  the  good  agreement  between 
the  calculated  and  tabulated  values  and  the  relative  ease  with 
which  a  desired  degree  of  accuracy  was  achieved,  the  approach 
using  approximating  series  was  employed. 

However,  as  indicated,  the  method 'had  to  be  modified. 

The  reason  for  this  is  evident  from  Figure'  3.  It  may  be 
seen  that  the  desired  probability  is  equal  to  the  area  under  the 
probability  curve  from  -  OO  to  t.  Because  of  the  symmetrical 
nature  of  the  normal  probability  curve,  the  following  technique 
was  applicable.  The  area  under  the  curve  was  obtained  for  the 
interval  -t  to  t.  This  value  was  halved  and  either  subtracted 
from  0.5  (for  negative  t)  cr  added  to  0.5  if  t  was  positive.  Thus, 
the  probability  for  -  CO  to  t  v/as  determined. 

In  addition,  the  subroutine  PROB  calculates  the  normal 
random  variable  associated  with  the  random  variable  x  from: 


where  A  is  the  mean  and  S  is  the  standard  deviation.  In  terms  of 
the  symbols  used  in  the  code: 


The  value  of  Z  is  adjusted  too,  using  the  expression: 


each  time  PROB  is  called.  This  is  necessary  since  Z.^  =  f(A,  S). 

The  corrections  for  S  and  A  are  calculated  using  the  current 
values  of  P,  Q,  and  Z  as  obtained  from  the  subroutine  PROB. 

The  new  values  of  S  and  A  are  used  to  evaluate  F^(S,  A)  and 
F£(S,  A}  and  then  compare  these  with  zero.  If  both  are  sufficiently 
close  to  zero,  the  order  of  accuracy  may  be  determined  by  the  User, 
For  instance,  1  x  10"^  and  1  x.  10“^  were  chosen  for  F^S,  A)  and 
F2(S,  A),  respectively,  for  the  code  shown  in  Appendix  A  while 
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THE  PROBABILITY  THAT  X  LIES  BETWEEN 
AND  u  IS  DEFINED'  AS: 

P(x)-mrf 

J -CO 


NORMAL  FREQUENCY  FUNCTION  FOR  THE 
NORMALIZED  RANDOM  VARIABLE  t 


FIGURE  3 


e 


1  x  10"^  and  l  x  10“^  were  chosen  for  the  code  of  Appendix  B. 
The  mean,  standard  deviation,  Fi(S,  A)  and  F^fS,  A)  are  then 
written  out  and  transfer  is  made  to  subroutine  SOLVE.  This 
subroutine  is  used  to  calculate  such  statistical  quantities  as  the 
standard  error  of  the  mean,  standard  error  of  the  standard 
deviation,  the  sure-firt  probability  point  and  a  corresponding 
safe  or  NO- FIRE  probability  level.  If  the  values  of  Fl  (S,  A)  and 
F2(S,  A)  are  not  sufficiently  close  to  zero,  however,  the 
iterations  are  continued  so  that  new  corrections  for  A  and  S  are 
obtained.  This  process  of  iteration  is  allowed  to  continue  until 
F  j(S,  A)  and  F£(S,  A)  become  as  small  as  desired  or  until  the 
number  of  iterations  reach  some  predetermined  value  chosen  by 
the  User  {1,000  iterations  were  chosen  in  the  code  shown  ir. 
Appendix  A,  100  for  the  code  shown  in  Appendix  B). 


The  values  of  Fj(S,  A),  F2(S,  A)  the  next  to  the  last  values  of 
S  and  A,  the  final  values  of  S  and  A  and  the  number  cf  iterations 
are  written  out  for  each  run.  This  is  done  so  that  if  all  1,000 
(or  100)  iterations  are  carried  out  because  the  error  in  either 
Fi(S,  A)  or  F2(S,  A)  is  larger  than  the  bound  specified  for  each, 
one  is  able  to  observe  two  successive  values  of  S  and  A  and 
decide  whether  convergence  has  actually  been  achieved.  The 
results  as  computed  by  each  code  for  three  tests  are  shown  at  the 
end  of  the  respective  code.  Note  that  for  these  three  tests  only 
two,  three  or  four  iterations  were  required.  This  rapid 
convergence  is  due  to  the  gc^d  approximation  of  the  mean  and 
standard  deviation  obtained  using  subroutine  FOREST, 


Once  the  maximum  likelihood  estimates  of  A  and  S  are 
approximated,  subroutine  SOLVEis  called  for  purposes  of 
estimating  the  standard  errors  associated  with  A  and  S  ar.d 
calculating  the  sure-fire  and  safe  or  NO- FIRE  energy  levels. 
To  estimate  the  standard  errors  two  parameters  G  and  H  must 
be  determined.  These  are  defined  by: 


and 


G  = 


Z.M /•£*.{ 

\  '  TV-/ 
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where: 


C  >  0 

0 


,  C^O  . 


Subroutine  SOLVE  makes  use  of  the  definitions  of  G  and  H 
and  calculates  these  two  quantities  directly.  Furthermore  it 
calculates  the  high  and  low  probability  points  and  the  other 
secondary  statistics  (standard  errors). 


The  standard  error  of  the  mean  is  calculated  from 
-  S6//  a7  ,  the  standard  error  of  the  standard  deviation 
from  r ^  ~  S  fi  /  'f/\J  *  the  standard  error  at  the  P 

percent  energy  level  from  %  *  G  *  *  A  SS  * 

In  the  above,  N  is  the  total  number  of  fires  ana  X  is  a  constant 
associated  with  a  particular  percent  point.  The  sure-fire  and 
safe  or  NO- FIRE  energy  levels  are  then  calculated  from: 


Sure-Fire  =  A  +  KS  (11) 

NO  FIRE  =  A  -  KS  (12) 


for  the  50%  confidence  band.  To  add  confidence  to  these  reliability 
points,  adjustment  to  Eauatmn  (11)  and  (2)  are  made  using  Sure- 
Fire  =  Sure-Fire  +  C  -.1  - s'and  NO- FIRE  =  NO- FIRE  -  a  N  +  C\t 


where  C  and  C‘  are  constants  associated  with  a  particular 
confidence  interval.  That  is,  the  estimates  of  the  sure  -fire  and 
NO- FIRE  energy  levels  should  be  correct  P  percent  of  the  time 
or  the  risk  of  being  wrong  is  anticipated  for  only  l-.P  percent  of 
the  time  (P  =  90,  95,  99, . , . ). 


As  stated  earlier,  the  code  in  Appendix  A  computes  all 
quantities  based  on  the  smaller  value  between  the  number  of  fires 
or  fails  obtained  in  the  Bruceton  test.  On  the  other  hand,  the 
code  in  Appendix  B  computes  all  quantities  without  regard  to  which 
value  is  smaller  except  in  weighing  averages  in  which  case  the 
smaller  is  more  heavily  favored. 

Finally,  the  main  program  writes  out  the  desired  statistical 
values 
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Regardless  of  which  code  is  used  the  data  cards  are  arranged 
identically.  The  following  describes  this  arrangement: 

Card  1  —  Six  fields  of  10  columns  each  are  reserved. 
Columns  1-10  and  11-20  are  reserved  for  constants  associated 
with  the  reliability  and  confidence  level,  respectively. 

Columns  21-30  and  31-40  are  reserved  for  specifying  the 
reliability  and  confidence  level,  respectively.  Columns  41-50 
are  reserved  for  specifying  the  NO- FIRE  energy  level.  Columns 
51-60  are  reserved  for  a  second  constant  associated  with  the 
confidence  level. 

Card  2  --  Four  fields,  one  of  20  columns,  one  of  5 
columns  and  2  of  10  columns  each  are  reserved.  The  first  field 
is  used  for  identification.  The  second  contains  the  number  of 
energy  levels  in  the  test.  The  third  and  fourth  contain  the 
lowest  energy  level  in  log  units  and  the  difference  between 
successive  energy  levels  (step  interval)  in  log  units,  respectively. 

Card  3  --  Twelve  fields,  arernating  between  eight  and 
four  columns  each  are  reserved.  The  fields  of  eight  columns 
each  contain  the  energy  levels  in  ascending  order.  The 
corresponding  number  of  fires  at  each  energy  level  is  specified 
in  the  right  adjacent  field  of  four  columns. 

Card  4  --  Six  fields  of  10  columns  each  are  reserved. 

The  number  of  fails  at  each  energy  level  --  starting  at  the  lowest 
energy  level  --  is  placed  in  these  fields  from  left  to  right. 

Table  II  on  Page  52  and  Page  20  and  21  of  Reference  4  contain 
the  constants  associated  with  several  reliability  and  confidence 
levels. 

This  example  will  demonstrate  the  input  data  cards  for  an 
actual  set  of  test  data.  The  table  below  shows  the  typical  results 
of  a  Bruceton  test.  It  is  desired  to  calculate  the  99.9%  sure-fire 
and  safe-fire  reliability  levels  at  the  95%  confidence  level: 
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Firing 

Current 

Log  of 

Firing  Current 

Test  Level 

Number  of 

Fires 

Number  of 
Fails 

374 

2.57287 

1 

0 

4 

392 

2.59329 

2 

4 

410 

2.61278 

3 

10 

10 

429 

2.63246 

4 

9 

2 

450 

2.65321 

5 

1 

0 

From  this  table  and  Reference  4  the  data  cards  for  input  may¬ 
be  prepared. 
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Card  1 

Card  Columns 
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The  output  using  these  data  cards  was  obtained  for  both  codes 
and  is  shown  at  the  end  of  each  (Appendix  A  and  B). 

The  code  of  Appendix  B  appears  to  be  superior.  First,  both 
the  number  of  fails  and  numbers  of  fires  are  utilized  in  computing 
the  various  statistical  quantities.  Thus,  a  larger  effective 
sample  size  is  employed  which  for  non-exact  measurements  is 
advantageous.  Also,  examination  of  the  results  at  the  end  of  each 
code  shows  that  a  higher  ALL-FIRE  and  lower  NO-FIRE  are 
obtained  using  Code  B  if  the  number  of  fires  and  fails  differ.  This 
is  desirable  for  increased  reliability  and  safey.  If  the  number  of 
fires  and  fails  are  equal,  the  results  using  either  code  are  identical. 
Consequently,  for  the  two  codes  use  of  that  of  Appendix  B  is 
advocated. 

In  certain  instances,  the  estimates  of  A  and  S  from  the  first 
two  moments  (approximation  method)  are  poor  approximations  to 
the  maximum  likelihood  values.  '  .  ,  In  particular  when 

the  step  size  between  test  levels  becomes  large  (d  >  2S)  and  when 
small  sample  sizes  are  used  (<50)  the  approximations  to  A  and  S 
become  poor  --  that  is,  the  maximum  likelihood  estimates  of  A 
and  S  are  not  well  approximated  by  the  first  two  moments. 

If  a  difference  in  the  two  estimates  does  exist,  the  maximum 
likelihood  quantities  should  be  used  because  they  are  expected  to 
be  better  than  those  generated  by  the  first  two  moments. 
Furthermore,  the  standard  errors  of  the  mean  and  standard 
deviation  are  based  on  their  respective  maximum  likelihood 
estimates  and  will  not  be  applicable  to  the  approximation  values 
when  a  large  difference  exists.  This  means  that  the  calculated 
extremes  (sure-fire  and  NO-FIRE  energy  levels)  will  be  poor 
estimates  of  the  true  quantities  as  well. 

To  determine  the  effect  of  sample  size  on  the  estimates  of 
the  standard  deviation,  a  series  of  tests  were  conducted  using 
sample  sizes  of  20,  30,  50  and  100  items.  The  mean  and  standard 
deviation  for  each  sample  were  calculated  by  estimating  the 
maximum  likelihood  quantities  and  by  utilizing  the  first  two 
moments.  Because  the  Bruceton  method  is  constructed  to  provide 
test  data  near  the  mean,  that  quantity  is  nearly  always  a  good 
estimate  of  the  population  value .  However,  the  standard  deviation 
is  not  so  well  determined.  Therefore,  in  what  follows,  special 
attention  is  paid  to  the  effect  of  sample  size  on  the  standard 
deviation  only. 
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Appendix  C  shows  calculated  values  of  the  mean  and  standard 
deviation  as  determined  by  the  maximum  likelihood  equations  and 
from  the  first  two  moments.  It  is  immediately  obvious  that  little 
difference  exists  between  the  means  as  calculated  by  either  method. 
But  for  small  samples  the  standard  deviation  as  calculated  from 
the  first  two  moments  differs  quite  appreciably  from  that  obtained 
utilizing  the  maximum  likelihood  equations.  As  stated,  standard 
errors  are  based  on  the  maximum  likelihood  quantities.  In  turn, 
the  extremes  of  the  probability  distribution  are  likewise  influenced 
because  of  their  dependence  on  the  mean,  standard  deviation  and 
associated  standard  errors. 

It  is  apparent  that  samples  containing  ever  100  items  should 
be  used.  The  use  of  maximum  likelihood  theory  for  the  cases  in 
which  the  sample  size  is  less  than  100  is  equally  advocated. 
Specifically,  Table  9  shows  that  for  samples  of  20  and  30  items 
large  differences  between  the  standard  deviations  occur  frequently. 
For  samples  of  size  50,  better  precision  is  achieved.  Neverthe¬ 
less,  some  poor  results  are  still  evident.  For  samples  of  size  100 
even  closer  agreement  between  the  two  methods  of  calculation  is 
achieved  and  it  is  anticipated  that  larger  sample  sizes  would 
produce  results  differing  only  slightly  from  one  method  to  the  other. 

It  should  be  pointed  out  that  the  good  results  in  the  50  item 
samples  were  due  primarily  to  a  knowledge  of  the  past  history  of 
the  test  vehicle.  However,  the  usual  situation  is  one  in  which 
little  or  no  data  is  available  for  the  item  since  it  is  the  purpose 
of  these  tests  to  generate  such  data.  In  these  'cases  it  is  imperative 
to  use  large  randomly  chosen  samples  or  if  this  is  not  practical, 
the  maximum  likelihood  equations  should  be  applied  in  calculating 
statistical  quantities. 

An  example  v/ill  best  illustrate  the  effect  of  using  small  sample 
data  in  calculating  the  extremes  of  a  population.  For  a  sample  of 
size  30,  these  were  the  results: 
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. . . . . — . . 

Maximum 

Likelihood  Equations 

Approximate 

Formulae 

Mean 

2.20095909 

2.18977 

Standard  Deviation 

0.18969359 

0.14878 

Standard  Error  of 
Mean 

0.041085 

0.03283 

Standard  Error  of 
Standard  Deviation 

0.073707 

0.05342 

Sure-Fire  Level 

3.27356932 

2.78843 

NO- FIRE  Level 

1. 12834880 

1.59111 

These  values  are  recorded  in  log  volts.  If  the  sure-fire 
and  NO- FIRE  levels  are  converted  to  volts,  the  potential  danger 
becomes  obvious  --  that  is,  using  maximum  likelihood  equations 
the  sure-fire  level  is  about  1,877  volts  as  against  615  volts  using 
the  approximation,  and  the  corresponding  NO-FIRE  levels  are 
13.5  volts  and  39  volts,  the  former  value  was  calculated  from 
the  maximum  likelihood  equations.  Consequently,  within  the 
probability  limits  indicated,  the  sure-fire  and  NO- FIRE  levels 
differ  roughly  by  a  factor  of  3,  Practically,  this  means  that  for 
this  case  in  which  a  firing  voltage  of  about  1,877  volts  is  required 
to  provide  sure-fire  conditions,  only  615  volts  might  be  made 
available  by  a  User  relying  on  approximate  formulas  only.  For 
similar  reasons,  safe  voltage  levels  may  be  inaccurately  specified. 
Thus,  anticipated  reliabilities  are  not  closely  achieved. 

Although  the  discrepancies  between  the  two  methods  become 
less  pronounced  with  increasing  sample  size,  maximum  likelihood 
methods  should  be  employed  in  view  of  their  assumed  superiority 
and  the  availability  of  a  digital  computer  algorithm. 
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APPENDICES 


APPENDIX  A 

FORTRAN  CODE  FOR  DETERMINING  THE  MEAN  AND  STANDARD  DEVIATION 
FROM  THE  SOLUTION  OF  THE  MAXIMUM  LIKELIHOOD  EQUATIONS  BY 
MEANS  OF  THE  METHOD  OF  SUCCESSIVE  APPROXIMATIONS 

(METHOD  A) 


BLANK  PAGE 
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BRUCE 


EFN 


SOURCE  STATEMENT 


IFN(S) 


C  BRUCE  STATIS  MAX  LIKE  ESTIMATORS 

DIMENSION  H(99),  Z (99 ) ,P ( 99) ,Q (99 ) , XN { 99 ) ,X( 99 ) , Y< 99 ) , W(99  ) 

DIMENSION,  DZSI.99),  DZA ( 99 ) , DOS ( 99 ) , DQA ( 99 ) ,DPS ( 99 ) , DPA( 99 ) 
DIMENSION  XNFI99)  ,ANTH(.50>  ,  I0ENT(2d) 

COMMON  GSQRD  ,  HSQRD  ,  SIGSQA  SIGSQS  ,  SIGMAP  ,  FIR50 
COMMON  FN0T5u»F IREK»FNOTK 

i*no.  of  Energy  levels  ,  s=step  sizeidist.  between  2  energy  levels 

A=*M IDDLE  ENERGY  LEVEL! IF  I  IS  ODD)  OR  AVG.  OF  MIDDLE  2  ENERGY  LEVE 
IIS  EVEN) 

I  READ  (5  ,2 )  RELK,CONKiRELV,C0NLV,XRELV,CST 

^•FORMAT {  6F10.0  ) 

CALL  FOREST  { XN , XNF , I , ANTH , «F I  RE . AFAI L , SFI RE, IDENT , SFA I L  ) 

DO  7855  J  *  1,1 
7855  H(J>  «  ALOGli(ANTHU)  ) 

H(J)*ENERGY  LEVELS  ,  XN ( J ) =N0.  OF  FIRES  AT  A  GIVEN  H(J)  ,  YlJ)»NOR 
1RAND0M  VARIABLE  ,  XU)*A8S.  VALUE  OF  Y(J)  ,  P(J)*PROB.  OF  FIRE  » 
QU»*PR0B.  OP  MO  FIRE  t  Z  (_4.i  *DERI  VATI  VE  0FPU» 

SUMX  *  0*0 
SUMNF  »  0.0 
DO  1216  J  *  1,1 
SUMN  *  SUMN  *XNU) 

1216  SUMNF  »  SUMNF. +XNFU) 

IF  (SUMN  -  SUMNF  )  8571,8571,1758 
8571  A  *  AFIRE 
S  *  SFIRE 
GO  TO  857-8 
1758  A  «  AFAIL 
S  »  SFA IL 

8-578  IF  (  SUMN.GT  .  SUMNF -r  CA..L  SWAP  (XN,XNF,X*) 

'24  WRITE  ( 6, 100  )  ( IDENT  U ) , J  =  1 ,20  ) 

100  FORMAT  (  1H1/1X,  2-jA2/30X  ,  66HCALC  •  OF  MEAN  AND  STD.  DEVIATION  BY  SOLN 
1.  OF  MAX.  LIKELIHOOD  EQNS.//) 

WRITE  (6  ,lu)’ 

10  FORMAT { 4QX , 12HENERGY  LEVEL , 25X , 18HN0.  OF  FIRES/FAILS//) 

WRITE  (6,11  HHU),XN(J),J  =  l,I) 

11  FORMAT  (  43X,F*10. 5,30X,  F4.0//  ) 

WRITE  (6  , 10 1 J ) A , S  ,  I 

1010  FORMAT ( 5 jX , 7HAVG.  »  , F 1 J. 5// 50X , 17HAPPRX  STND  DEV  *  , F 10. 5//50X , 23 
1HN0.  OF  ENERGY  LEVELS  *  ,14//  ) 

COUNT  OF  THE  NO.  OF  ITERATIONS 
ITCNT  »  0 

6  CALL  PROB  (H, X,  Y,  Z,  A ,  S  ,  >», Q  ,  I ) 

BEGIN  SUCCESSIVE  APPROXIMATION SINEW TON  RAPHSON  METHOD) 

DERIVATIVE  OF  Z  WITH  RESPECT  TO  S  AND  A 

DERIVATIVE  OF  Q  WITH  RESPECT  TO  S  'A"ND  A 

DERIVATIVE  OF  P  WITH  RESPECT  T0  S  AND  A 

DO  436  J»li I 

DZS( J)*t (H( J)-A)**2-S*S)/S»*3*Z( J) 

DZA(J  )»(H(J)-A) /!S*S)*Z  J  ) 

DQS (J  )  =-( H!J )-A )/( Q IJ ) *  ») *  *2*Z ( J ) 

DWA(J)»-Z(J)/(g(J)«*2*S) 

DPSU)a(H(J)-A)/(P(J)»S)»*2*Z(J) 

OPA(J)=ZU)/(P(J)*P(J)*S) 

436  CONTINUE 
SUMlaO.O 
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EFN  SfURCE  STATEMENT 


IFNCS ) 


SUK2=0.0 
SUM 3* 0.0 
SUM4*C.O 
SUM5*C.G 
SUM6*0.0 
SUM *0.0 
SSUM*0.0 
00  555  J*2,I^* 

-SUHl*SUMl*xAj|JJ*(0ZS(  J-lI/t  (  J-l  J  *005(3-1  )*Z(J-I)-DZS(Ji /Pi  JI-CPSCJ 
1)*Z(J)J  A 

SUM2»SUM2*Xfcgjl*(Z(J-n/Q(  J-ll-ZIJJ/PlJ)) 

SUM3«SUM3*xM J I • ( OZ A ( J-l ? /c  C  J-i 1 ♦OQA  ( J-l >  *Z  C J-l l-OZAl JI/PI J >-0P A( J 
1)*Z(JI) 

SUM4*SUM4+XA^J)*(  (H(  J-l)-A)*OZSt  J-l ) /C( J-l )*IH(  J-l  l-AI  »CQSC  J-l » •ZC. 
1J-1)-(H( JJ-AJ*DZS( J)/P( J)-(H(J|-A)*DPS£ J)*Z(J») 

SUH5*SUH5*X#ofj i * { (HI J-1I-A)*Z( J-l » /QU-1 1- IH( J)-A)»Z( JJ/P( J#  J 
SUM6=SUM6*XN(|  J ) • ( H ( J-l ) *OZA ( J-l J /O I J-l » ♦« ( J-l I *OQA C  J-l J  «Z ( J-l ) -Z ( J 
1-ll/QI J-ll-A/QI J-1)*0ZA( J-I  )-A*I(  J-1!«DQA(  J-ll-Hi JI/PI J»«DZA( JI-DP 
2A(  JI*H(  J]*Z(  J)*Z(JI/P(J)*A*Z(  J)*DPA(  J)*A/«>(  J)*DZA(J)) 
SUM*SUM+XWJJ«(Z(  J-il/Q(J-l)-Z(  JI/PI  Jl) 

SSUM*3SUMVj^jj|*({H( J-1I-A)*Z( J-l)/OCJ-l)-(H(JI-AI*ZCJ)/P(Jil 
555  CONTINUE 

C  DERIVATIVE  OF  F1(S,A)  WITH  RESPECT  TOS  AND  A 

C  DERIVATIVE  3F  F2(S»A)  WITH  RESPECT  TO  l  AN 0  A 

C  FAS=Fi(S,A)  ,  SAS*F2(S,AI 

OFOS*1.0/S*SUH1-1.D/S**2  *SUH2 
0FDA*1.0/S*SUM3 

0S0S*1.C/S**2  *SUM4-2.0/S**3  *SUH5 

0S0A=1.0/S**2  *SUM6 
FAS=»1  • J/S*SUM 
SAS*l.u/S**2*SSUK 
303  D*0F0S*DSDA-DFDA*0S0S 
C  CORRECTION  FOR  S 

XH*(OFOA*SAS-DSOA*FAS)/P 
C  CORRECTION  FOR  A 

XK*( DSDS*FAS-DFOS*SAS i /o 
C  CORRECTED  A 

AA«A*XK 

C  CORRECTED  S 

SS*S+XH 
17  S“SS 
A*AA 

C  VALUE  OF  F 1  {  A ,  S )  AND  F2(A,S)  OBTAINED  FROM  FINAL  A  AND  S 
ZER0i«O.O 
ZEK02*O .0 
00  69  J  =  i:  *  I 

ZER01«ZERO1*XNC  J  J  *  (  Z  C  J-l) /(.I  J-  1)-Z(JI/PU1I 
69  ZER02*ZER02>XN< J)*( (H( J-l)-A)«  Z ( J-l ) /Q i J-i )-(H( J )-A ) »Z ( J ) /P ( J )) 
IF(AOS<ZER0l)-l.E-3J  68,68,15 
68  IF(ABS(ZER02)-l.E-4)  16,16,15 

15  ITCNT  *  ITCNT  ♦  1  • 

I F (  ITCNT  -  1000  )  151,16,16 
151  G0  T 0  6 

1 6  I  T  E  (6  ,6666)ZER01,Zf  R02 , A , AA , S , SS 

6666  F0RMATI SOX, 8HZER01  »  ,  E20 .  V.  lib  OX ,  8HZER02  *  ,  C20 . 8// 50X  ,  9HNLME  AN  « 
1,F10.8//50X, 7HMEAN  *  , F10. 8//50X , 1 5HNLSTND.  DEV.  *  ,F10.8//50X, 
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BRUCE  -  EFN  S0UK.CE  STATEMENT  -  IFNCSI  - 

221HSTAN0ARD  DEVIA7I3N  » 

WRITE  I 6  » ill ) I TCNT 

111  F0R«ATt5OX,2OH.N3.  0F  ITERATirNS  =  ,15//) 

CALL  S0LVECP ,Q,XN,Y,Z,i,A,S,H.RELK,C3NK,CST) 

WRITE  (6  ,77) 

77  F0RMATI4DX,49HCALCULATI0N  0F  ALL-FIRE  AND  N8—FIRE  ENERGY  LEVELS//) 
WRITE  (6  ,971  iGSQRD,HSQRO,  SIGSCJA,  S IGSQS  »SIGMAP,  FI  R5D*  FN3T50 

l,R£LV,C0NLV,*FIREX,XRELV,C3NLV,rf.?TK 
971  F0RHAT I 53X.11HG SQUARED  *  , FI  *.  6//5';X  ,  1  jHHSQUARED  «,F10.6//50X,29HS 
1TAND.  ERRBR  3F  MEAN  SQRD.  *  , F  l'j.6// 53X,  36HSTAND.  ERRBfl  0F  STAND. 
2DEV.  SCRO.  *  ,F 1Q.6//5CX, 35HSTAND.  EKR3R  AT  P  PERCENT  HEIGHT  »  ,F1 
30.8//5*»X, 33H99.9FIRE  AT  5DC3NFIDENCE  LEVEL  *  ,F10.8//50X,32HQ. 1FIR 
4£  AT  5QC3NFIOENCE  LEVEL  *  ,F10. 8//50X.F6.2 ,7HFIR£  AT.F6.2, 19HC0NFI 
50ENCE  LEVEL  »  ,F1D.8//5DX,F6.2,7HFIRE  AT.F6.2, 19HC0NF1DENCE  LEVEL 
6>  , F10.8//  ) 

G0  T0  1 
END 
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8R0P 


EFN  S0URCE  STATEMENT 


I FN ( S  > 


SUBR3UTINE  PR081  H,X, Y,Z,A,S,P,Q, i) 

OIKENSI0N  •  H{99)»X{99),Y{99),Z(99)tP199)»Q(99) 

00  7  J  *  1*  I 
HMASQ  »  ( H( J )— A ) *«2 
EXPNT  *  HHASU/l 2.»S#S) 

Z(J)  *  ( 1 • /2. 50663) /EXP (EXPNT) 

Y(J)  «  (HU  I— A)  /S 
X(J)  *  ABSCY(J) ) 

IF  (XU  1-1*20  )  21,21,22 

21  XSO  *  X(J!*X(J) 

P(J) *0.79788455*XU)*(0. 999997 74-XSQ»(0.16659433-XSQ»(  .24638310E-0 
il-XSQ*.23974867E-02) ) )  - 
IF  (  Y  C  J  )  )  400,  401,  402 
403  P  (  J  )  *  0.5-  -  P  (  J  )  /  2.0 
G0  T0  433 

401  P  (  J  )  *  0.5 
G0  TO  433 

402  P( J)*?U)/2.0*0.5 

403  Q  (J  )  *  1.3-  P  (  J  ) 

G0  T0  7 

22  IF  ( X ( J )  -  2.90  I  25,26,26 

25  XSQ  *  X(J)  *  XU) 

PU)  «  i.O 

PT  *1.0 
FACT  *  1.0 
BOO  *  3.0 

29  PT  *  -PT*XSQ/(2.0*FACT  ) 

T  *  PT/000 
P(J)  *  P(J)  ♦  T 

IF  (A3SIT)  -  0.33037  )  63  ,70  ,70 
70  FACT  *  FACT  ♦  1.3 
000  *  000  ♦  2.3 
G0  T0  29 

63  P(J)  *  0.79788455*XU)*PU) 

IF  (  Y  ,(  J)  )  500,  501,  502 
503  P  i  J  )  =  0.5  -  P-  {  J  )  /  2.0 
G0  T0  503 

501  P  (  J  )  =  0.5 
G0  T0  503 

502  P (J ) *PU ) /2.0+3.5 

503  Q  1J  )  «  1.0-  P  {  J  ) 

G0  T0  7 

26  RXSQ  *  1.0/IXU  )  *XU )  ) 

P( J)  *  1.0- 3.7978 C455*EXP(-XU)*XU)/2.3)/XU)*(  1 .0-RXSO* ( 1 .0**RXSQ 
1*  ( 3.0-RXSQM  15.0-RXS01C5.0)) ) ) 

IF  (  Y  (  J  )  )  600,  601,  602 
603  P  (  J  )  *  0.5  -  P  {  J  )  /  2.0 
00  T0  603 

601  P  (  J  )  =  0.5 
G0  10  603 

602  P(  J  )  *P  (  J>/2. 0+0.5 

603  0  ( J  )  *  1.0-  P  {  J  } 

7  CONTINUE 

RETURN 

ENO 
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F0RTES 


£}  N  SOURCE  STATEMENT 


IFN(S) 


SUBR0UT INE  F0RESTlFiX,FAX»N,X,UfIX,UFAX»VFIX»IDENT,VFAX  )  L 
ALCG(W)=.4343*a'L0G(W) 

ANT(W)=EXP(VW/.4343) 

01  MENS  1 0NX  (  50)  ,  1D£NT(20)  , F  i  X < 99 ) , FA*  199)  ,S  (30,15) 

Data  S  /o4» , 25. ,40. ,43. ,46. ,49. , 51. , 16*0. ,2 

1 5.  *27.  ,31. .,35.  ,37.  ,42.  ,4'  ,49. ,  5  1 .  ,  12*0 .  ♦  25 . ,  26  .  , 28  .  ,  29 . ,  31 .  ,  3 

24. . 36. .39. .40. .43. .46. .48. .4 9. .31. .5.0. ,25. ,27., 20. ,31., 32. ,34., 36 

3», 37., 39., 49 . « 41 • , 42 . *  44 • ,46. ,43. ,45. ,51. ,4*0. ,25. ,26. ,27., 2/., 28. 
4, 3j . , 3j . ,31. ,32. ,33. ,35- ,  36. ,37. ,34. ,41. ,42. ,43. ,44. ,46., 49 . , 49  .  , 5 

51. . 25. .26. .27. .27. .28. .29. .39. .30. .31. .32. .33. .34. .35. .36.. 37. .38. 

6 . 40 .. 41.. 4 2. .43. .45.. 46. ,47^,49., 51., 52. ,25. ,25., 26., 27., 29., 29., 3 

70.  ,30. ,  32.  ,32.  ,33., 34.  ,:>5.  ,36.  ,  38.  ,39.  ,40.  ,42.  ,43.  ,44.  ,45.  ,46. ,  47. 
.8,49., 51., 52. ,25., 26. ,27. ,27. ,29. ,29. ,30. ,31., 32., 33. ,34.%  35., 36.,  3 

97. .  38.. 39.. 40. .41. .43. .44. .45. .47.  ,49.  ,50.  ,51. ,53. ,26. ,26. ,27. , 28. 
A, 29., 30., 31., 31., 32. ,34. ,3 5., 36.  ,37.  , 39. , 39 . , 40 . , 4 1  . , 42 . , 44 . , 44  . , 4 

65. . 47. .49. .50. .52. .53.. 26. .27. .20. .28. .29.. 31. .31. .33. .34. .34. .36. 
C ,  37.  ,38.  ,39.  ,49., 41., 42. ,43.,  44  .  ,45.  ,46..,  47  .,  49  .,51. ,52. ,53. ,26. ,2 
07., 28., 29., 30., 31. ,32. ,33. ,34.  ,35. ,36. ,37.  , 38. , 39. , 40 . , 41 . , 42 . , 43 . 
E,44.,45. ,46. ,47. ,49. ,50. ,52. ,54./ 

READ  (5,639) (1DENT(J),J-1,20) ,N,C,D 
_639  F0RMAT (  20A1,  1 5, 2F  10.0)  .. 

READ( 5, 637 )  ( X { I ) , F I  X (I) , I *1 ,N ) 

_637  F0KMAT  (  6  (  F8  •  0 ,  F4.0  ).).. 

REA0(  5,633)  ( FAX  C I) , I *1 , N ) 

..638  .F0RMAT (  6F10.0  ) 

AF I X*0. 0 
AFAX  =  0.0 
8FIX*0.0 
.BF AX*0 • 0 
SF  l  X*0 • 0 
SFAX«0.0 
1*0 
3010  W= I 

AFIX=AFIX*W*FIX( 1+1  ) 

AFAX*AFAX+W*FAX( 1+1) 

BF IX  *  BFIX  +  ( W**2  )  *  PIXU  +  l) 

8FAX*BFAX+(U**2)*FAXU  +  1) 

SF I X  =  SF  I  X^F I  X ( l+l) 

SFAX=SFAX+FAX( Ui) 

1  =  1  +  1 

IF ( I-N+1J3010, 3010, 1115 
1115  R  *  D 

TFIX*( (SFIX*BrIX)-AFIX**2)/(SFIX**2) 

TFAX* { ( SFAX*BFAX)-AFAX**2) /(SFAX«*2) 

13  UFIX*C+R*(AFlX/SFIX-0.5) 

UFAX=C+K*( AFAX/SFAX*0.5) 

I F { TF I X-0 • 3 )  118,23,23 
118  C  I  X=TF 1 X*100 • 

L I X=C  I X 

4111  WEEIX  «  500. 

1=0 

112  WEEIX  »  AMINKWEEIX,  A8S(UF  IX-ALCG(X(  I  +  1  i  )  )  ) 

1  =  1  +  1 

IFU-N+1)  112,112,113  35 

113  RIX*A8S(WEEIX/RI 


F0RTES  -  EFN  S3URCE  STATEMENT 

KIX  *  I F I X  ( R I X  *  22.  )  ♦  1 
VFiX*{  S ( L I X»  K l X ) / 1 00 • )*R 
2424  IF( TFAX-O. 3) 11? 1 2525 1 2525 
117  CAX*TFAX»iuO • 

LAX*CAX 
114  WEEAX*500. 

I  *  j 

115.  WEEAX  *  AMIN1(W££AX,ABS<UFAX-ALCG{X(  I  +  i>)  )  ) 
1*1+1 

I F ( I-N+ 1)115*115*116 
116  RAX*ABS(WE£AX/R) 

KAX  *  I F IX  { RAX*22. ) +1 
VFAX  *.(  S ( LAX s KAX) /10G. ) *R 
G2  TV  26 

23  VFIX=1.62*R*(TFIX+C. 32 4 i « SORT l $r l X/JSFIX-l. 
G0  T0  2424' 

2S25  VFAX* 1, 62*R* ( TFAX+C . 329 ) *  SORT ( SFAX/ { SFAX-1 . 


RETURN 

FND 


IFN(S) 
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f 

PAWS  —  EFN  SOURCE  STATEMENT  -  1PM  (S) 

Subroutine  SWAP  (XN,  XNF,  1) 

Dimension  XN(99),  XNP(99) 

D0  4143  J  -  2,  I 
4143  XN(J)  -  XNF  (J-l) 

RETURN 

END 
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VALVE 


SOURCE  STATEMENT 


1FN(S) 


-  EFN 


SUBROUTINE  SOLVE  {P,Q,X.\i,Y,Z,  I  ,A,S.HtRELK,C0NK,CST  ) 

C0MM0N  GSQRD t H-SQRDt  S IGXQA «  SIGSQS * S IGMAP, F  I R50  »  FN0T5Ot  FIREK ,  FN0TK 
OI^ENSJL3fiL  H{99),Pt99),U99)  ,XN(99) , Y(99)  ,Z(99)tW(99) 

SUMNJ  =  0.0 
Will  *  1.0 
WJSUM  *  J.O 
GSUM  *  0.0 
HSUM  »  0.0 
00  51  J*2» I 

51  W(J»aQU-l)/p(J-l)*W(  J-l) 

00  61  J*l,l  __ 

SUMNJ  *  SUMNJ  ♦  XN(  J ) 

61  WJSUM  >  WJSUM  ♦  W(J) 

00  741  J*2#I 

GSUM  «  GSUM  ♦  W<  J-i  )  *  (  Z  ( J-l  )**2/Q(  J  *1 )  **2«-Z  ( J )  **2/P(  J )  **2  ) 

T-rl  HSUM  a  HSUM+WI J-l ) *  I Y{ J-l )  **2*ZI J-l  I**2/QI J~1  )**2*Y  ( J ) **2*Z ( J )**2/ 
1P(JJ**2) 

GSQRO  * WJSUM/GSUM 

HSQRO  *  WJSUM/HSUM 

SIGSQA*S*S/SUMNJ*GSQRD 

SIGSQS  *  S*S/SUMNj*HSOPa 

SIGMAP  *  SQRT(.SIGSGA+KfcLK**2*$IG$QS  ) 

FIR6J  *r  A-*-RELK»S 
FN0T5O "a  A-RELK*S 

FIREK  *  F IR5O+C0NK* ( { SUMN J+CST ) /SUMNJ ) *SI GMAP 
FNfTK  =  FN0T5O-C0NK»{ ( SUMNJ+CSTJ /SUMNJ ) *S I GMAP 
RETURN 
'END 
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M2  Squib  Test  D 


Calc*  bf  mean  and  sto.  oeviatibn 

ENERGY-  LEVEL 
2.57287 

2.59329 

2.61278 

2.63246 


BY  S0LN.  0F  MAX.  LIKELIH00O  EONS. 

N0.  0F  FIRES/FAILS 

O. 

4. 

10. 

9. 


2.65321  i, 

AVGi  *  2.60870 

APPRX  STNO  OEV  *  0.02159 

N0.  0F  ENERGY  LEVELS  »  5 

2ER01  *  -0.52526593E-04 

2ER02  *  0.28066803E-04 

NLMEAN  »  2.60874474 
MEAN  *  2.60874474 
NISTND.  06V.  *  0.02079978 
STANOARD  OEV I  AT  I 0N  »  0.02079978 
Ntf •  0F  ITERATIONS  *  2 

CALCULATION  0F  ALL-FIRE  AND  N0-FIRE  ENERGY  LEVELS 
GSQUARED  «  0.696013 

HSQUAR60  «  0.791929 

STAND.  ERROR  0F  MEAN  SQRD.  »  0.00001254 
STANO.  ERR0R  0F  STAND.  DEV.  SORD.  *  0.00001428 
STAND.  ERROR  AT  P  PERCENT  HEIGHT  «  0.01220044 
99.9FIRE  AT  5OC0NF l OENCE  LEVEL  «  2.67301604 
0. IF  IRE  AT  5OC0NFIOENCE  LEVEL  »  2.54447341 
99.90FIRE  AT  95  .OOC0NF I DENCE  LEVEL  *-  2.69812453 
O.IOFIRE  AT  95.OOC0NF I DENCE  LEVEL  »  2.51936489 
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MS  Can  Test  A 


CALC.  0F  MEAN  AND  STD. 

ENERGY  LEVEL 
0.14301 

0. 16435 

0.18469 

0.20683 


0F  MAX.  L1KELIH00D  EQNS 
N0.  0F  FIRES/FAILS 
C. 

2. 

9. 

10. 


DEVIATION  BY  S0LN 


0.22531  3* 

AVG.  *  0.18602 

APPRX  STNO  DEV  *  0.02279 

N0.  0F  ENERGY  LEVELS  *  '  5 

ZER01  *  0.19371510E-05 

ZER02  »  0.22351742E-07 

NLMEAN  *  0.18685731 
MEAN  »  0.18685731 
NLSfNO.  DEV.  »  0.02369443 
STANDARD  DEV  I AT  1 0N  *  0. 0236£443‘ 

•M0.  0F  ITERATIONS  «  3 

CALCUL ' TI 0N  0F  ALL-FIRE  AND  N0-FIRE  ENERGY  LEVELS 
GSQUARED  *  0.846374 

HSQUARED  »  1.498057 

STAND.  ERROR  0F  MEAN  SQRD.  *'  0.00001980 
STANO.  ERR0R  0F  STAND.  DEV.  SQRD.  «  0.00003504 
STAND.  ERR0R  AT  P  PERCENT  HEIGHT  *  0.01882551 
99.9FIRE  AT  5OC0NF I DENCE  LEVEL  »  0.26007312 
j. IFIRE  AT  6OC0NF I OENCE  LEVEL  *  0.11364149 
99.90FIRE  AT  95.00C0NF I DENCE  LEVEL  *  0.29881602 
0.10FIRE  AT  95 .OGC0NF I DENCE  LEVEL  *  0.07489859 
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M6  Cap  Test  E 


$»- 


CALC . 


0F  MEAN  AND  STD.  DEVIATION  BY  S0LN. 
ENERGY  LEVEL 
2.40993 


0F  MAX.  L1KELIH00O  EONS. 
N0.  0F  FIRES/FAILS 
0. 


2.42488 


2. 


2.44091 


9. 


2.45484  12. 

2.46982  2. 

AVG.  =  2.44130 

APPRX  ST NO  OEV  =  0.01477 

0F  ENERGY  LEVELS  «  5 

Z.ER01  =  0 .2676  248  6E-04 

ZER02  =  0 .751 11166E-06 

NLMEAN  =  2.4*135469 

mean  =  2.44135469 

NLSTNO.  OEV.  =  0.01453041 
STANDARD  DEVIATION  *  0.01453041 
NO.  0F  I TERATI 0NS  =  2 

CALCULATION  0F  ALL-FIRE  AND  N0-FIRE  ENERGY  LEVELS 
GSQUARED  =  0.785916 

HSQUAREO  »  1.171202 

STAND.  ERROR  0F  MEAN  SQRD.  *  0.00000663 
STAND.  ERROR  0F  STAND.  DEV.  SQRD.  »  0.00000989 
STAND.  ERROR  AT  P  PERCENT  HEIGHT  *  0.01005381 
99.9FIRE  AT  5OC0NF I DENCE  LEVEL  *  2.48625368 
0. IFIRE  AT  5  3C0NF I OENCE  LEVEL  *  2.39645568 
99.90FIRE  AT  95 .QOC0NF I DENCE  LEVEL  »  2.50690502 
0.10FIRE  AT  95 .OOC0NF I OENCE  LEVEL  *  2.37580431 
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APPENDIX  B 


FORTRAN  CODE  FOR  DETERMINING  THE  MEAN  AND  STANDARD  DEVIATION 
FROM  THE  SOLUTION  OF  THE  MAXIMUM  LIKELIHOOD  EQUATIONS  BY 
MEANS  OF  THE  METHOD  OF  SUCCESSIVE  APPROXIMATIONS 

(METHOD  B) 


oooo  o  ooo  ooo 


BRUCE 


£FN  L0URC6  STATEMENT 


I CN( S ) 


C  BRUCE  STATIS  MAX  LIKE  EST I MAT0RS 

DIMENSION  H ( 99 ) »  Z {9^ ) , P { 99) ,0(99 ) » XN( 99 ) • X { 99 ) « Y( 99 ) , W ( 99 »  . 

DIMENSION  *DZS  ( 99 )  ,  DZA {  99  ) , DQS ( 99  )  , D0A( 99  ) , DPS( 99  )  , DPA( 99  ) 

D I  MENS 1 0N  XNF ( 99) » ANTH { 50 ) » IDENT  1 20 ) 

COMMON  GSQRD  ,  HSQRD  ,  SIGSQA  ,  SIGSQS.j.  SIGMAP  FIR50 
C0MM0N  FN0T5O.FIREK, FN0TK 

I=N0..  0F  ENERGY  LEVELS  ,  S*STEP  S  I ZE  ( DI  ST. _  BETWEEN  2  ENERGY  LEVELS. 
A=MIDDLE  ENERGY  lEVELIIF  I  IS  3DD)  0R  AVG.  0F  MIDDLE  2  ENERGY  LEVE 
IIS  EVEN) 

1  READ  (5  .2)  RELK , C0NK ,R£LV ,C0NLV , XRELV ,CST 

.2  F0RMATI  6F10.0  ) 

CALL  FOREST(XN,XNF,I,ANTH,A,S,IDENT  ) 

00  7855  J  =  1*1 
7855  H ( J )  =  AL0G1O(ANTH( J) ) 

H(J)=ENERGY  LEVELS  ,  XN(J)-N0.  0F  FIRES  AT  A  GIVEN  H(J)  »  Y(J)'N0P. 
1RAND0M  VARIABLE  ,  X(J)*ABS.  VALUE  0F  Y(J)  ,  P(J)«PR0B.  0F  FIRE  » 
Q(J)=PR0B.  0F  N0  FIRE  ,  Z ( J ) =DER I  VAT  I VE  0F  P(J) 

24  WRITE  (6,100  ) ( IDENT ( J  ) , J~  1 ,20  ) 

100  FORMAT ( 1H1/1X,20A2/30X, 66!  iCALC .  0F  MEAN  AND  STD.,  DEVIATION  BY  S0LN 
1.  0F  MAX.  LIKELIH00D  EQNS. // ) 

WRITE  (6  ,10) 

13  FORMAT ( 40X, 12HENERGY  LEVEL ,25X, 18HN0.  0F  FIRES/FAILS//) 

WRITE  (6,11  )  (HtJ),XN(J)  ,  XNF  ( J)  ,  J  *  1,1  ). 

11  F0RMAT (40X,FlG.5»33X,F4.Qt3X»F4.0//) 

WRITE  (6  , 1 010 ) A , S  ,  I 

1310  F0RMAT(5OX,7HAVG.  *  ,F 10. 5//>UX , 1 7HAPPRX  STND  DEV  *  » F 10. 5//50X , 23 
1HN0.  0F  ENERGY  LEVELS  *  ,14//  ) 

C0UNT  0F  THE  N0 .  0F  ITERATIONS 
__  I TCNT  =  0 

6  CALL  PR08  (H , X, Y, Z , A , S , P,U ,  I ) 

BEGIN  SUCCESSIVE  APPR0X IMATI 0NS ( NE WTON  RAPHS0N  ME7H0D) 

DERIVATIVE  0F  Z  WITH  RESPECT  T0  S  AND  A 

DERIVATIVE  0F  Q  WITH  RF.SPb.CT  T0  S  AND  A 

DERIVATIVE  0F  P  WITH  RESPECT  T0  S  AND  A 

D0  436  J*l, I 

DZS(J )  =  ( (H(J)-A)**2-S*S)/S*«  i*Z ( J  » 

OZA(J)*(H( J)-A)/(S*S)*Z( J) 

DQS(J)=-(H(Ji~A)/(Q(J)*S)**2*Z(J) 

DQA( J)«-Z(J)/(Q(J)**2*S) 

DPS(J)a(H(J)-A)/(P-(J)*S)**2»Z(J) 

OPA( J)»Z( J)/(P( J)*P( J)*S) 

436  CONTINUE 
SUM1=0.G 
SUM2*0.0 
SUM3=0.0 
SUM4*0 . 0 
SUM5=0 . 0 
SUM6  =  0  *  0 
SUM=0 . 0 
Sr>UM  =  0  •  0 
D0  555  J  =* 2  , 1 

SUM1*SUMU-XNF(  J-l)MDZS(J-l)/Q(J-iUOQSU-l)*Z(  J-l)  )-XM  J )  •  (  DZS  ( J  ) 
i/P ( J ) +DPS ( J ) *  Z ( J  )  ) 

SUM2«SUM2  +  XNF(J-1)*(Z( J-l)/Q( J-l) ) -XN ( J ) * ( Z ( J ) / P ( J  )  ) 

SUM3»SUM3+XNF( J-l > * ( DZA { J - 1 > /0( J~1 ) +  (>QA( J-l ) »Z ( J-l  ) )-XN( J)* ( DZA( J ) 
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l'/PC  JI+DPAi  JI*Zl  JJ) 

SUM4*SUK4*XNF( J-l)*(  IHI J-1)-4J*0ZS( J-ll/GI J-1)*IH( J-l)-AI«OQSI J-1J 
1*Z( J-l)  )-XN( J )•  I (H{  J)-A)»DZi I Jl/Pl  J)+IK(  J>-A)«DPSC  J)*Z  C  J)  * 
SUM5=SUM5*XNFCJ-1>*I  1H(J-.1)-A>»Z(J-1)/QI  J-m-XN(J)«{(HfJ!-A)«Z(Jl 
I/PIJ)) 

SUM6»SUM6*X?iF ( J-l ) • I H ( J-l ) «(-ZA ( J-l ) /Q( J-l > ♦«{  J-l  J *DQA  I  J-l  I *Z ' J-l l-_ 
1ZCJ-l)/Q(J-l)-A/aiJ-l)*DZA<J-i)-A.ZIJ-l)*0QAiJ-l>)-XNlJ)*'HlJ|/PtJ 
2 ) *OZA C  J ) +OPA l J ) *H I J ) *  Z I J ) - Z ! J ) /? ( J ) -A*Z I J  >  *OPA ( J ) -A/P ( J ) *0Z A ( J  J  )  . 
SUK»SUM*XNf* J-ll*IZi J-1I/9I J-1))-XN(J)*IZC J)/P( Jil 

SSUM*SSUM*XNF(J-i)M  CHCJ-1)-A)«ZC  J-l ) /OI J-l M-XNC Jl *i f HI J)-A)«Z! J ) 
i/?i  jn 
555  CONTINUE 

DERIVATIVE  OF  FliS.A)  WITH  RESPECT  TBS  AND  A 
DERlVAl 1VE  BF  F2IS.41  WITH  RESPECT  T0  S  AND  A 
FAS*F I I S, A )  ,  SAS*F2(S.A) 

0FDS*l#0/S«SUMl-1.0/S**2  *SUM2 
OFDA»l. j/S«SUH3 

DS0S*1.0/S»*2  *SUM4-2.G/S**3  *SUM5 

DSDA*1.0/S**2  -sim 
FAS*1.0/S*SUM 
SAS*1.0/S»*2»SSUM 
303  D=DFPS»OSDA-DFDA*DSDS 
CORRECTION  F0R  S 
XH»(DFDA*SAS*-OSDA*FAS!/D 
CORRECTION  F0R  A 
XK*( DSOS«FAS-OFDS*SAS ) /D 
CORRECTED  A 
AA*A>XK 
CORRECTED  S 
SS*S+XH 
IT  S=SS 
A*AA 

VALUE  OF  F1(A,$)  AND  F2IA,S)  OBTAINED  FROM  FINAL  A  AND  S 
ZERO l =0.0 
ZER02O.O 
00  69  J*2 »  I 

ZER  .=ZER01*XNF(J-1)»(Z(J-I)/U(J-1 ) )-XN( J ) * IZI J l/Pl Jl > 

69  ZER02»ZER02*XNF(J-l)*»  )-AI*Z( J-l )/Q( J-l)  l-XNt  II*  I  IH(J)-A)*ZI 

IJI/PI J) ) 

IFIABSIZEROIl-l .E-4)  68,68,15 
68  IF(ABS(ZER02)-l.E-5)  .6,16,15 

15  ITCNT  =  ITCNT  ♦  1 

I F ( ITCNT  -  100  )  151,16,16 

151  G0  TO  6 

16  WRITE  16  ,6666)ZER01,ZER02,A,AA,S,SS 

5666  FORMATl  5DX,8H76R01  *  , E20. 8//5'  X , 8H?C R02  =  ,E20.8//50X ,9HNL*EAN  * 
l » F i J. 8/ /50X » .7HMCAN  =  ,Fl0.8//5._  X,1  5H.NLS TND •  DEV.  *  ,F10.8//50X, 

22 1HST ANOARD  DEVIATION  *  ,F10.8//> 

WRITE  16  ,111) ITCNT 

111.  FORMAT!  SOX,  20.1N0*  OF  ITERATIONS  =  ,15//) 

CALL  S0LVEIP,O,XN,Y,Z, I ,A,S,H,RELK,C0NK,CST,XNF  ) 

WRITE  (6  ,77) 

77  c0RMAH  4OX.,A9HCALCULATI0N  OF  ALL-FIRE  AND  N0-FIRE  -ENERGY  LEVELS//) 
WkITE  16  ,971)GSQRD,H3&RD,SlGSQA,SIGSQStSIGMAP.FIR50,  FN0T50 

i,RELV,C0NLV,FIREK,XRELV,C0NLV,FNOTK 
971. FORMAT* 50X  1  lHGSQ'JARED  *  ,  F 1 0. 6//5DX  ,  lOHI  SQUARED  =  ,  F10 . 6//50X  ,  29HS 


ts 


0RUCE  -  £FN  S3URCE  STATEMENT  -  IFNI SI  - 

HAND*  ERR0R  0F  MEAN  SURO.  *  ,F 10. ~ //50X,36HSTAND.  ERR3R  0F  STANC. 
20EV.  SQRO.  *  »Fi0.S//50X,25HSTAN0.  ERROR  AT  P  PERCENT  HEIGHT  *  ,F1 
30 •  8// 50X « 33H99 •  9F*I  RE  AT  5C  CO'lF  IDENCE  LEVEL  *  ,F10.8//53X,  32h0.  IFIR 
A£  AT  50C8NFIDENCE  LEVEL  *  ,Fi ;.8//50X,F6.2,7HFIRE  AT, F6.2, 19HC0NFI 
50ENCE  LEVEL  *  • F10.8//50X, F6.2 ,?HF IRE  AT,F6.2»19HC0NFIO£NCE  LEVEL 
6*  , F10.8//  » 

G0  T0  1 
END 


45 


BK0P 


EFN  S0U1CE  STATEMENT 


IFNIS) 


Subroutine  probi  h,x,y,z,a,s,p,c,i  j 

DIMENSION  ’  HI99),XI99),YI99I,Z199),P599),Q199) 

00  7  J  *  1,  I  • 

HMASQ  *  IH(JI-AJ**2 
EXPNT  *  HM ASQ/ 1 2  •  *  S*  S ) 

ZlJ)  *  ( l./2.50663)/EXPIEXPJ*T) 

YIJJ  *  (HUJ-AJ/S 
XIJI  *  ABSI Y I J  i ) 

IF  IXIJ)-1.20  )  21,21,22 

21  XSQ  *  X(J)*X(J) 

PI  J)  =0. 79783453*XI J  >  *  ( 0.99999774-XS0*  CO.  16659433-XSQ*  (.•246.38310Er.0_ 
ll-XS0*.23974367E-02)  ) ) 

IF  I  Y  I  J  )  )  409,  401,  402 

400  P  I  J  I  *  0.5  -  P  I  J  :  *  2.3 
G0  TO  403 

401  P  I  J  )  *  0.5 
G0  T0  433 

402  PI J)=Pl J J/2.G+0.5 

403  Q  (J  )  <  1.3-  P  I  J  ) 

G0  T0  7 

22  IF  IXU)  -  2.90  I  25,26,26 

25  XSQ  *  XIJ)  •  XIJ) 

PIJ)  *  1.9 

PT  =  1.0 
FACT  *  1.0 
0DD  *  3.0 

29  PT  *  -PT*XSQ/ l 2 .0* FACT  ) 

• T  =  PT/0OD 
PIJ)  =  PIJ)  ♦  T 

IF  ( APS ( T )  -  9.00G37  )  63  ,70  ,70 
10  FACT  *  FACT  ♦  1.0 
BOD  =  000  ♦  2.0 
G0  T0  29 

63  PIJ)  *  0.79788455*XI J)*P( J) 

IF  l  Y  l  J)  )  500,  501,  502 

500  P(J)=0.5-PIJ)/  2.3 
G0  T0  503 

501  P  l  J  )  =  0.5 
G0  T0  503 

502  PI  J)=P( J)/2.0+0.5 

503  Q  (J  I  *  1,5-  PIJ) 

G0  T0  7 

26  RXSQ  =  1.0/IXI J)*XI J) ) 

PIJ)  =  1.0-3.79788455*EXPI-X{J)*X{J)/2.)>/X{J)*{ 1 .O-RXSQ* I 1 .0-KXSQ 
i»  { 3. O-RXSQ* I15.G-RXSQ*1C5.0) ) ) ) 

IF  I  Y  I  J  )  )  600,  601,  602 

600  PIJ)*  0.5  -PIJ)/  2.J 
G0  T 0  603 

601  P  I  J  )  =  0.5 
G0  T 0  603 

602  PSJ)=PlJ)/2.0+0.5 

603  Q  I J  )  *  i .0-  P  I  J  ) 

7  CONTINUE 

RETURN 

END 


F0RTES 
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IFN(S) 


SUBROUTINE  F0R£ST(FIX,FAX,N,X,U,V,IDENT  ) 

ALCG(U) =.4343*AL0G ( W ) 

ANT(Wl*EXP(W/.4343) 

OIMENS 10NX ( 50)  , I DENT (20),F!X(99) , FAX ( 99) , S  (30,15) 

DATA  S  /64*0., 25. >40. ,43. ,46. ,49. ,51. ,16*0. ,2 


15. ,27., 

31 

..,35 

.,37., 42. ,44 

.  .47 

.,49 

.,51. ,12*0 

.,25 

•  • 

26  . , 

28., 

29 

•  9 

31 

..3 

24., 

36., 

39 

.  ,  40 

. , 43. , 46 . , 48 

.,49 

.  .51 

•  1 9*0 • 1 25 ♦ 

,27. 

,28. ,31., 32 . 

♦  34. 

,36 

3., 37. ,39. 

,  4"). 

, 41 • , 42 . , 44. 

,  46. 

,  48  • 

,49., 51., 4*0 . , 

25 

.,26 

.,27 

•  f 

27 

•  9 

28. 

4,30 

.,30 

•  9 

31., 

32. ,33. ,35., 

36., 

37., 

39. ,41. ,42 

.,43 

•  V 

44., 

46  .  , 

49 

*  , 

49 

.,5 

51., 

25., 

26 

.,27 

.,27.  ,28.  ,.’9 

.  ,30 

.,30 

.,31. ,32., 

33.. 

34 

.,35 

.,36 

•  t 

37 

•  9 

38. 

6*40 

•» ,  41 

•  t 

42. ,43. ,45* ,46. , 

47., 

49., 

51. ,52. ,25 

.,25 

•  f 

26  . , 

27., 

29 

.*  , 

29 

.,3 

70., 

30., 

32 

.,32 

.,33. ,34. ,35 

.  ,36 

.  ,  38 

.,39., 40 . , 

42., 

43 

..44 

.,45 

•  t 

46 

♦  9 

47. 

_8, 49 

.,51 

•  t 

52., 

25. ,26., 27., 

27., 

29., 

25. ,30. ,31 

.,32 

•  , 

33., 

34., 

35 

•  , 

36 

.,3 

97., 

38., 

39 

.,40 

.,41. ,43., 44 

.,45 

.,47 

•  , 49 . , 50 .  , 

51., 

53 

.  ,  26 

.  ,  26 

•  t 

27 

•  9 

28. 

A, 29 

.,30 

•  t 

31., 

31. ,32. ,34., 

35., 

36., 

37. ,39. ,39 

.,40 

•  , 

41., 

42., 

44 

• , 

44 

.  ,  4 

B5.  » 

47., 

49 

.,50 

.,52. ,53. ,26 

.,27 

.,28 

•  , 28 . ,29.  , 

31., 

31 

.,33 

.,34 

•  t 

34 

•  9 

36. 

C,  37 

.  ,  38 

•  9 

39., 

40. ,41. ,42., 

43., 

44., 

45.  ,46.., 47 

.*49 

.  , 

51., 

52., 

53 

•  , 

26 

• ,  2 

07., 

28., 

29 

.,30 

.,31. ,32. ,33 

.,34 

.,35 

.  , 36 . , 37 . , 

38., 

39 

.  ,  40 

.,41 

•  9 

42 

•  9 

43. 

E»44 

.,45 

•  t 

46.  , 

47  .  ,  49  .  ,  5G_«  , 

52., 

54./ 

READ  (5.639)  (IDENT( J ) , J=i, 20)  ,N,C,D 
639  F0RMAT( 20A1.I5.2F1 0.0) 

READ( 5, 637 )  ( X { I ) , F I  X ( I ) . 1 =1 .N) 

63T  ,F^RMAT(6(F8.0,F4.0) ) 

READ! 5. 638 )  ( FAX ( I ) , I *1 , N) 

638  F0KMAT(  6F1Q.0  ) 

AFIX»0.0 
AFAX  *0.0 
BFIX=0.0 
BFAX*0.0 
SFIX*0.0 
SFAX*0.0 
1*0 
3010  W*  t 

AFl  .  AFIX*W*FIX( I+l  ) 

,AFAX*AFA.X+W*FAX(  I  +  l ) 

8FIX  *  BFIX  ♦  ( W**2  )  *  FIX(I+1) 
BFAX*BFAX+(W**2)*FAX( I  +  l) 

SF I X*SF I X+F I X( I  +  l) 

SFAX*SFAX+FAX( I+l) 

1*1  +  1 

IF( I-N+l) 3010, 3010, 1115 
1115  R  *  D 

TF I X*  (.(  SFI  X*BF  I  X  )-AF  IX**2  )  /  ( SF  I  X**2 ) 

TFAX*( (SFAX*8FAX)-AFAX**2)/(SFAX**2) 

13  UFlX*C+R*(AFiX/SFlX-0.5) 
UFAX*C+R*(AFAX/$FAX+0.5) 

. IFITF lX-0.3)  118,23:23 
118  C IX*TF 1 X*100» 

1 1 X*C  I X 

4111  WEE l X  *  500. 

1*0 

112  WEEIX  *  AMIN1(WEEIX,ABS(0F IX-ALCG(X( I+l) )) j 
1*1  +  1 

IFU-N+D112, 112,113 

113  R I X*A8S { WEE ! X/R ) 
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SUBROUTINE  SOLVE  IP ,Q» XN, Y , Z , I , A , S »H , RELK, CBNK.CST.XNF  ) 

COMMON  GSQRD, HSQRD, SIGSQA,SIGSQS,3IGMAP,FIR50,  FN0T5O,FIREK,FN0TK 
DIMENSION  H{99),P(99)»QI99),XN(99)»Y(99>  , Z 199)  »H(99 )  *XNFf 991 
SUMNF J =3.0 
.SUMNJ  *  0.0 
WJSUM  *  O.'J 
GSUM  *  0.0 
HSUM  *  0.0 
IZERO  *  i/2  ♦  i 
IM1  »  IZERO  -  1 
00  623  J  *  1,1 
JM1  =  J  -  1 
W1J)  *  1.0 

IF(H{ J )  -  HI IZERO  I)  358.347,223 
358  00  6262  K  *  J,IM1 
6262  WCJ)  »  W  ( J )  ♦  PIK)  /  Q(K) 

GO  TO  623 
347  WIJ1  =  1.0 
GO  TO  623 

223  DO  727  L  *  IZERO  ,  JH1 
727  WIJ)  *  HU)  *  QIL)  /_  P I L) 

623  CONTINUE 

00  741  J  *  1,1 
WJSUM  *  WJSUM  ♦  W(J) 

SUMNF J=SUMNF J+XNF ( J ) 

741  SUMNJ  =  SUMNJ  ♦  XNIJ) 

DO  61  J  =  2,  I 

GSUM  =  GSUM  ♦  HI J-i)*CZ(J-l)**2/QI J-l )  **2«-Z  1 J  >  **2/P  I J )  **2  > 

61.  HSUM  =  HSUM+WI  J-i)*(Y(  J-1)**2*ZC  J-l)**2/Q(  J_-U**2t.Y.(J).**2*Z(  J)**2/ 
1P( J)**2) 

GSQRO  =  WJSUM/GSUM 

HSQRD  =  WJSUM/HSUM 

SIGSQA=S»S/SUMNJ*GSQRO 

SIGSOS  *  S*S/$UMNJ*HSQRD 

SGSQAF=S*S/$UKNFJ*GSQRQ 

SGSQSF=S*S/SUMNFJ* HSQRD 

SIGMAP  =  SORT ( SIGSQA+RELK**2*S IGSQS  ) 

SGMAPFsSQRTI SGSQAF+RELK*«2«  SGSQSF ) 

.FIR50  *  A+RELX*S 
FN0T5O  =  A-RELK*S 

_F  IREK  =  FIR5O+C0NK*((SUMNJ  +  CST)/SUMNJ)*SIGMAP 
FN0TK  =  FN0T5O-C0NK* ( ( SUMNJ+CST) /SUMNJ )*SIGMAP 
_FIREKF  =  FIR5O+C0NK»( ( SUMNF  J  +  C ST ) /SUMNF J) *SGMAPF 
FN0TKF  =  FN0T5O~C0NK* { ( SUMNF J+CST) /SUMNF J ) *SGMAPF 
S I  GSQA=SQRT  (  (  S  I  GSQA**2  *-SGSQAF»  *2  )  /?. .  ) 

S I  GSQS  =  SQRT  (  (  SI GSQS*»2  *-SGSQSF*»2)/2.  ) 

F ! REK  =  SQR  T { { F I REK»*2+F I REKF »*2 ) /2 . ) 

FN0TK  3  SORT ( FNUTK  *  FNJTKF  ) 

RETURN 
END 
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KIX  =  IFIX  (RIX  *  22.  *  ♦  1 
VF1X»(  S( L IX t KIX) /ICO. ) *R 
2424  lF(TFAX-0.3)llV, 2525, 2525 
117  CAX*TFAX*lOO. 

LAX*CAX 

.114  W£EAX*530. 

I*u 

115  WEEAX  =  AMIN1 { W£EAX, ABS { UF AX-ALCG ( X JI ♦! ) ) ) ) 

1*1  +  1 

IF( I-N+1J 115,115,116 

116  RAX*ABS { HEEAX/R ) 

KAX  *  IFIX  ( RAX*22. ) +1 
VFAX  *.<  St  LAX, KAX) /IOC. ) *R 
G0  T0  26 

23  VFIX*1.62*R*(TF IX+0.C29 )»SORT (SFIX/(SFIX-1.)) 

G0  T0  2424 

2525  VFAX=1. 62*R* (TFAX+O. 029)* SORT < SFAX/ ( SFAX-1 . ) ) 

26  U*( tSFAX*UFAX)  +  (SFIX*UFIX) ) / ( SFAX  +  SF I X ) 

V*SQRTUSFAX*(VFAX**2)+SFIX*(VFIX**2) )/(SFAX+SFlX) ) 

RETURN 

END 


M2  Squib  Test  D 


CALC.  IF  MEAN  ANO  StD.  DEVIATION  BY  S0LN.  i.F  MAX.  LIKELIH00O  EQNS. 

. 

ENERGY  LEVEL  Nfl.  FIRES/FAILS 

2,57287  °*  *• 

2.59329  *•  *<>- 


2.61278 


10.  10. 


2,63266  9-  '2. 

2.65321  1-  0. 

AVG.  *  2.60967 

APPRX  STNO  OEV  »  0.02289 

N0.  3F  ENERGY  LEVELS  *  5 

Z.ERC  1 .  *  C,  1.6 19942  5.E-06 

ZER02  ■  .0.4401 66 3.3.E-0 6 

NLME  AN  _  2  . 610942.16 
MEAN  »  2.61094216 
N.LSTNO,..  OEV  ...  *  ..0_.02_41510.3_ 

STANDARD. DEVIATION  *  0.02415103 
N0,  0F  I  I'i_R  AT  I  0NS  *  3 

CALCULATION  0F  ALL-FIRE  ANO  N0-F.IR.E  ENERGY  .LEVELS 
GSOUAREO  «  0.8 263.8.3 

HSOUAREO  -  _1.496004 

ST ANO^  ERROR  0F  MEAN  SQRO.  *  0.00001932 
*TAND.  ERR0R.  0F  STAND.  OEV.  SQRO.  *  0.00003498 
STANO.  ERR0R  AT  ?  PERCENT  HEIGHT  ■■  0^01916319 
99.9FIRE  AT  5OC0NFIOENCE  LEVEL  •  2.68556884 
1. _1F IRE  AT  SOC0NF.I&ENCE  LEVEL  *  2.53631544 
99.90FIRE  AT  95 .OOC0NF IOENCF  LEVEL  *  2.72416380 
O.iOFIRE  AT  95 .OOC0NF  t  OENCE  LEVEL  «  2.49772042 


50 


MS  Cap  Test  A 


CALC.  ,0F  MEAN  AND  STD.  DEVIATION 
ENERGY  LrVEL 
G.143'31 

0.16415 

0.184<>9 

0.20603 


BY  S0LN.  0F  MAX.  LIKELIH00D  EQNS. 

N0.  0F  FIRES/FAILS 

0.  2. 

3.  9. 

10.  10. 

11.  3. 


0.22531  3.  0. 

AVG.  «  0.18533 

APPRX  STNO  OEV  «  0.02337 

N0.  0F  ENERGY  LEVELS  •  5 

ZER01  *  0.15720725E-05  . 

ZER02  »  0.10360964E-07 

NLM6AN  =  0.18430348 
MEAN  «  0.18430348 
NLSTNO.  OEV.  =  0.02588342 
STANOARO  DEVIATI0N  ■  0.02588342 
N0*  0F  ITERATI 0NS  *  4 

CALCULATION  0F  ALL-FIRE  ANO  N0-FIRE  ENERGY  LEVELS 
GS0UAR60  =  0.835367 

HSQUAREO  *  1.510502 

STAND.  ERROR  0F  MEAN  SQRO.  «  0.00002206 
STAND.  ERROR  0F  STAND.  OEV.  SQRD.  0.00003989 
STANO.  ERROR  AT  P  PERCENT  HEIGHT  *  0.01945744 
99.9FIRE  AT  5OC0NFIDENCE  LEVEL  «  0.26423325 
0-.1FIRE  AT  5OC0NF I OENCE  LEVEL  *  0.10432371 
99.90FIRE  AT  95 .OOC0NF I OENCE  LEVEL  «  0.305433 
O.IOFIRE  AT  95 .OOC0NF l OENCE  LEVEL  «  0.063157 
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.  M6  C«p  leatJS 


2.46982  2«  0. 

_AVG.  *  2.44130. 

APPRX  STNO  DEV  *  .0.01477 

_N0.  0F  ENERGY  LEVELS  *  5 

2ER01  *  0.26669353E-04 

ZER02  -  0.75314892E-06 

NLMEAN  ,*  2.44135469 
MEAN  *  2.4413.5469 
NLSTNO.  OEV._  *  „0. 01.4530.41 
^STANDARD  .DEVIATION  ■  0.01453041 
N0.  0F  ITERATIONS  *  2 

CALCULATI  ON  0F  .ALL-FIRE  AND  N0-FIRE  ENERGY  LEVELS 
_GS.QUA.REO.*  .0,7859 16. 

HSQUARED  1.171202 

.STAND.  .ERROR  0F_  .MEAN. SQRO..  •  0.00000663 
STAND.  ERR0R  0F  STAND.  DEV.  SQRO.  *  0.00000989 
STAND.  ERR0R  AT  P  PERCENT  HEIGHT  *  0.01005,381 
,_99.9FIRE  AT  5OC0NF I OENCE  LEVEL  *  2k48625368 
.0. IFIRE  AT  5OC0NF IOENCE  LEVEL  -»  2.39645568 
99.9DFIRE  AT  95 .OOCfNP I DENCE  LEVEL  *  2.50690502 
. 0 . 10F IRE  AT  95 .OOCINF IDENCE  LEVEL  *  2.37580431 
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